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Abstract 

The general quasiparticle propagator in dense quark matter is de- 
rived for equal mass quarks. Specialized to an NJL model, this propa- 
gator includes one new condensate, A3, in addition to the usual CFL 
condensate, Ai. The gap equation is solved in two NJL models and 
the dependence on the quark mass of the condensates and the gap is 
presented. Analytic approximations for the condensates are obtained 
and compared to exact numerical solutions. The results are shown 
to differ from those obtained by neglecting A3, especially for smaller 
values of A\. The two different NJL models presented are also shown 
to give different results when A3 is not neglected. The methods used 
in this paper can be generalized to the physical case where only the 
strange quark is significantly massive. 

I Introduction 



The physics of strongly interacting matter at high densities and low temper- 
atures has been the subject of much research in recent years. It has long been 
known [Jl], that at sufficiently high densities a system of quarks should form a 
condensate of Cooper pairs which breaks the SUc{<$) symmetry and becomes 
a color superconductor. The formation of the diquark condensates leads to 
gaps in the quasi-particle spectrum. The authors of |], @] estimated that the 
gaps were of the order of A ~ 10 _3 /i where /1 is the quark chemical potential. 
More recently it was shown at realistic values of \x in an instanton induced 
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NJL[] model that gaps of the order of 10-100 MeV could be obtained^. This 
stimulated a great deal of research^ in the ensuing years and has resulted 
in a proliferation of predicted superconducting ground states. These states 
may be realized in the cores of neutron or quark stars and lead to observable 
effects!, i I- 



It is widely accepted[10, 11] that for three colors and flavors the color 



superconducting ground state at asymptotic densities is the Color Flavor 



Locked (CFL) state|10 



tfpHfi) = a33(^-4o + a 66 (^ + 4^) (i) 

\A=2,5,7 

where the Greek indices are color indices, the Latin indices are flavor indices, 
the 3 and 6 subscripts refer to anti-triplet or sextet configurations in color 
and flavor spaces respectively and X A are the Gell-Mann matrices. At lower 
densities (/i ~ m^/4A), it is likely that the ground state is a superconducting 
state involving the condensation of Cooper pairs in the u — d sector only 
(2SC): 

(qlC 1 5 q^ = A 2SC e aP3 e^, (2) 

Finally at still lower densities the favored ground state will be ordinary 
hadronic matter. Two new phases have recently been predicted: Crys- 
talline Color Superconductivity |L2| and CFL with meson condensation fi"3|, 
\A . These predictions, while not necessarily at odds with one another, indi- 
cate that the transition region between the CFL state and hadronic matter 
is not completely understood. 

A standard way of studying color superconductivity is to solve the gap 
equation in the Nambu Gorkov formalism which is reviewed in Appendix A. 
The inclusion of a strange quark mass in the gap equation introduces two 
sets of complications: 1) massive quarks means that there are 4 new allowed 
Dirac structures for the condensates, and coupling between the condensates; 
2) the fact that the strange quark is different from the other quarks means 
that condensates involving the strange quark should be different from those 
with zero strangeness. 

1 The designation NJL (Nambu- Jona-Lasinio) model is used in this paper to describe 
models where the quark interaction is taken to be a four-fermion interaction at a point. 
2 For extensive references see the review articles ||, |j| . 
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In order to understand the implications of these two complications it is 
useful to separate them and understand them individually before tackling the 
full problem. In a previous paper by this author |0| the second problem of 
non-degenerate quarks was studied in the case where all quarks are massless 
but the strange quark is given a different chemical potential than the other 
two quarks. This paper concentrates on the first complication by considering 
the problem of equal mass quarks. 

Color superconductivity with massive quarks has been studied before. 
The papers |16|-||23|| have dealt with inclusion of quark mass or strange quark 
mass in NJL models. The papers |I3| , [H], ^3]| have dealt with inclusion of 
quark mass or strange quark mass in effective Lagrangian models. Finally [25. 
PB[ have dealt with this problem in the perturbative approach. All of these 
important papers have dealt with the problem at some level of approximation. 
The new results in this paper are solutions for all allowed diquark condensates 
applicable over a wide range of quark masses and coupling constants. The 
quasiparticle propagator derived herein is the general form for the case of 
equal mass quarks. It is shown that for the NJL model there is one and only 
one additional diquark condensate, A 3 , and it is included in the calculations. 
The approximate analytic solutions presented in this paper make numerical 
solution of gap equations for different parameter values unnecessary. Finally 
two different NJL models are analyzed in exactly the same framework so that 
they can be compared and contrasted. 

The authors of ]TB|| did discuss the fact that A 3 , which they refer to as 
k, must be non-zero but since it is small they neglect it in their calculations. 
The authors of include a similar diquark condensate in the 2SC phase, 
but their analysis does not indicate that no other diquark condensates would 
arise. In it was shown that A 3 can be large in an NJL model using a 
Lorentz non-invariant interaction and therefore have a significant effect on 
the gap. In this paper it is shown that, even for small values of A3, neglecting 
this condensate in the calculation can have a non-trivial effect on the gap. 

The results of this paper are an approximation to the Nf = 3 CFL case 
as the A 66 gaps are neglected to simplify the presentation. The results are 
exact for Nf = 2 and for the 2SC phase with Nf = 3. Therefore there is 
only one color-flavor structure and we drop the subscripts "33" and "2SC" 
as they are unnecessary. 

The general quasiparticle propagator is presented for the case of equal 
mass quarks. This propagator is then specialized to the NJL model where it 
is shown that there is only one new diquark condensate, A 3 , in addition to 
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the CFL condensate, Aj. 

The gap equation was solved in two different NJL models. The first model 
used is the simplest possible NJL model equivalent to scalar exchange. The 
second model used is an NJL model with the color structure of single gluon 
exchange. The second model is perhaps more physically motivated and is 
the model usually used in NJL analyses of color superconductivity. However, 
there is no reason to exclude scalar interaction terms and restrict the analysis 
to interactions motivated by single gluon exchange as a general low energy 
Lagrangian could include both. This being said it is interesting to see how 
much of an effect the specific NJL model has on the results. 

The gap equation was solved in the scalar NJL model and results are 
presented for the general set of condensates as a function of the quark mass. 
Approximate analytic solutions for the condensates are obtained and shown 
to agree well with exact numerical solutions. The results are compared with 
results obtained neglecting the effect of the new condensate A3. It is shown 
that the inclusion of A 3 alters the results especially for lower values of A^. 
The gap in the quasiparticle spectrum is shown to be a function of both A x 
and A 3 and inclusion of A 3 is even more relevant. 

The same analysis was performed in an NJL model with the color flavor 
structure of single gluon exchange. The results in this model are also altered 
by the inclusion of A 3 especially for lower values of A x but they have the 
opposite effect as in the scalar NJL model. 

These results are relevant to the analysis of where the gap equa- 
tion is solved neglecting A3. Their results are shown to be a reasonable 
approximation for the parameters they chose, but can be generalized using 
the techniques of this paper. The effects of chiral symmetry breaking have 



not been taken into account in this paper as in [^, ^TJ but such an extension 
would not be difficult. 

These results illustrate the effect of using the general form of the quasi- 
particle propagator in the solution of the gap equations for massive but de- 
generate quarks. The results are most significant for smaller values of the 
CFL condensate and they depend strongly on the specific NJL model used. 

The methods used in this paper are a significant new result on their own. 
They can be generalized to determine the general quasiparticle propagator 
and the general set of gap equations for the physical case where the up and 
down quarks are essentially massless, but the strange quark is massive. As 
well, this analysis could be repeated using perturbation theory in order to 
compare with a model that is exact at higher densities and energies. 
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The outline of this paper is as follows. In section II a basis for the Dirac 
structure of the condensates is presented which facilitates the rest of the 
analysis. In section III the general quasiparticle propagator is derived and 
then specialized to the NJL model. In section IV the gap equations are solved 
in the scalar NJL model. In section V the gap equations are solved in an 
NJL model which has the color structure of single gluon exchange. Section VI 
concludes with a summary of the results. Appendix A contains a brief review 
of the Nambu-Gorkov formalism. Appendix B gives the decomposition from 
@ of the gap matrix in terms of the products of Dirac matrices. 



II Basis for the Condensates 

The CFL condensate in ([I]) has the Dirac structure 75 and is the only spin 
zero condensate for the massless case in an NJL model. The most general set 
of spin zero condensates is given in Appendix B and includes eight different 
possible structures. In an NJL model the structures involving 7 • k automat- 
ically vanish leaving four possible Dirac structures: 75, 7570, 7o, I- Solving 
the gap equations give the result that there is only one new condensate, A 3 
with the Dirac structure: 7570- 

The analysis could be carried out in a basis of Dirac structures, but the 
equations would be very complicated. A crucial step in finding the general 
quasiparticle propagator and gap equations is choosing a basis which simpli- 
fies the analysis enough to make it tractable. 

In the massless case an arbitrary condensate matrix can be decomposed 



using a basis of orthogonal projection operators [gTfl: 

PZ(k)=A e (k)P h (k), e,h = ±l, (3) 
which are products of positive and negative energy projectors: 

A .(S) = 1±£|±« e = ±l, (4) 
and positive and negative helicity projectors: 

r 1 + e7 5 7o7- k 
p h{k) = h = ±1. (5) 

In this basis the condensate matrix is: 

a = A t n- (6) 

e,/i,±l 
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Similarly all objects, such as the bare quark and the quasiparticle propagator 
can be written in this basis. The basis of orthogonal projection matrices is 
extremely useful since it reduces products of these objects to the simplest 
possible form. 

In the case of massive quarks things are more complicated. The energy 
projector can be generalized to|T 



1 + e ( /?7o7 ■ k + a7o ) 
A e (A:) = y — 2 ^ e = ±l, (7) 

where (3 = \k\/Ek and a = mj E^. A set of analogous operators in the 
massive case given in [27| is: 

P c e h (k) = P c A e (k)P h (k), e,h = ±l c = r,l, (8) 

where: 

c = r,J, (9) 

is the chirality projector which projects onto right and left handed spinors. 
The definition r = + and / = — in analogy with the other projectors will 
be used in what follows. These operators are really quasiprojectors with the 
general product rule: 

P! h {k)P e J h> {k) = 5 eel 5 cc ,5 hh ,P e ch {k) - l -ece'c'{l - e'ch (3) S hh , P^(k). (10) 

They are not projectors because chirality and energy projectors do not com- 
mute: 

[P c ,A e ] = -eca 7 o75. (11) 
This basis is complete and the general gap matrix can be written[27|: 



A= E ^ h pe c h - (12) 

e,c,/i,±l 

This basis still involves many complications. 

The calculation is simplified by defining a new basis of true projectors 
and nilpotent operators: 

pm= E p° h (k) = p h (k) A%k), (is) 

c=-l,l 
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Qt(k)= E (ech-(3)P: h (k) = (eh 1 5 -(3)Pak), 

c=-l,l 



that satisfy the following relations: 

P e h {k)P^{k) - 
Qt(k)P^(k) - 
P e h {k)Qi{k) - 

Qt(k)QUk) 



See'S, 
5 



= —a 



ee' 

ee'$hh'Qh e (k), 
2 S-ee'8hh>Ph £ - 



(14) 



(15) 
(16) 
(17) 
(18) 



The Ph{k) are orthogonal projectors and the Q%{k) are nilpotent even though 
products of P%(k) and Q%(k) or Q e h (k) and Q e h \ mix the operators. P%(k) is 
an hermitian operator and (Q e h (k)Y = —Q~ h e (k). In the chiral limit (a — > 0, 
/3 — > 1) P%(k) are identical to the massless projectors and Q%(k) — > 0. Notice 
that the operators with different helicities are completely decoupled. 
The multiplication table for each helicity is: 





Pit 


Qt 


Ph 


Ql 


Pit 


(P h + 








Ql \ 


Qt 


Qt 








-c?Ph 


P h 





QX 


Ph 





Ql 


V o 


~o?P h + 


Ql 


o J 



(19) 



This set of eight operators is the set of operators on the space of solu- 
tions of the Dirac equation which gives the most sparse multiplication table. 
One would prefer a basis of only orthogonal projection operators as in the 
massless case, but such a set does not exist. Dirac spinors are 4 compo- 
nent complex vectors and as such live in 1Z 8 . One can definitely find eight 
orthogonal projectors in this space, but Dirac spinors must be solutions of 
the Dirac equation which places constraints on the components and reduces 
the dimensionality of the space. Therefore, you cannot have 8 orthogonal 
projectors onto the space of solutions of the Dirac equation. On the other 
hand you need 8 operators to form a complete set with which to construct 
all 8 independent Dirac structures. In the massless case, the left and right 
handed spinors decoupled so that each was a two component complex spinor 
living in TZ 4 with two constraints from the Dirac equation. It is therefore 
possible to find a set of 4 orthogonal projectors on this space. 
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The operators given above have the following representation as 2 x 2 
matrices: 



«-( -a!!).«-(!! o)- <»> 

This supports the reasoning above that there are only 4 independent pro- 
jectors since the operators for each handedness operate on a separate two 
dimensional space. This representation will be useful in what follows. 
This set of operators is complete as shown by the following relations: 



I = 


e,h=-l,l 


(22) 


7 5 7 ° 7 . k = 


e,h=-l,l 


(23) 


a 7 • k = 


- E QhW, 

e,h=—l,l 


(24) 


a 7 5 7° = 


E hQ e h (k), 

e,h=-l,l 


(25) 


7°7 • k = 


Y e(Q e h (k)+f3P^k)), 

e,h=-l,l 

E eh(QUk)+(3P^(k)), 

e,h=-l,l 


(26) 


7 5 = 


(27) 


a 7 ° = 


Y e(a 2 pe(k)-[3Ql(k)), 

e,h=-l,l 


(28) 


a 7 5 7 • k = 


- Y eh(a 2 PZ(k)-PQt(k)). 

e,h=-l,l 


(29) 



It should be noted that these equations are still valid in the chiral limit (albeit 
half of them are trivial). 

The gap matrix is decomposed in terms of these operators as: 

A= E (pa + c?j>t)PZ(k) + (? h -PP h )Q° h (k), (30) 

e,h=-l,l 

At = E (pa+» 2 ^)pz(k)-(e h -(3r h )Qr(k). (si) 

e,h=—l,l 



All objects in this analysis can be represented using this basis. The sparse 
multiplication table simplifies much of the analysis in the rest of this paper. 
In the next section the general quasiparticle propagator is derived using this 
basis. 



Ill General Quasiparticle Propagator 

The quasiparticle propagator is determined by the equation [|27f p| 

^{[G^-A^A^ 1 . (32) 

where A + = A, A~ = 7 A + 7 and Gq is the bare charge-conjugate particle 
propagator: 



G (k) = (fkv-rf-M) 



-i _ — ^7° + M 
~ {ko _ ^2 _ El 

(k -fi)- £ fc /? 7 °7 • k + E k aj°) 



(33) 



(k - /i) 2 - El 
o (jg -ti)+E k Ee.^-i.1 e [(a 2 - (3 2 )P e h (k) - 2(3 Q e h (k)} 
(ko ~ /i) 2 - El 



Defining: 



n 



G + (k) ((A; - / i)7 - 7 .fc + M)~ 
£ (AlP^k) + B e h Ql(k)) = 



e,h=-l,l 

the coefficients are: 

At = (kof - e e (a,r h y 

using the notation: 



h=-lA 



2ea 2 E k 



e\2 



(k - n) - eE k 
2eE k 



(34) 



(35) 



(k - fi) + eEi 



e ± (i^)= (E kT ^ 2 + e + oiU 2 



1/2 



(36) 



3 see Appendix A for a short review of the Nambu-Gorkov formalism. 
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n _1 = E Jo-Urw-Brt, (37) 



Using the 2x2 matrix representation of and Q e h given above: 

E - JL - 

; ,ht^i,i det fi/, 
where: 

detfi, = k*-2k*(El + ^ + (ti) 2 + (^) 2 + a^f + a\^) (38) 
+ 2E k k ((^+) 2 - te) 2 ) 

+ - /i) 2 (£ fe + /i) 2 + + /.) 2 (£+) 2 + (E k - /i) 2 (^) 2 

- a 2 (E k - pt)(E k + /,) ((^) 2 - (^,) 2 ) + + o?^l) 

and the definition: 



E k = Vl^P + m 2 , (39) 

has been used. The poles of the quasiparticle propagator are given by the 
roots of this quartic equation in k Q . The roots have well known analytic 
forms that are nonetheless quite complicated and are not necessary for the 
purposes of this paper. 

From (^4j) is it clear that: 

G + (k) = fi- 1 ((A; - /i) 7 ° - 7 • ^ + M) , (40) 

and therefore the general quasiparticle propagator is: 

G + (k) = (41) 
a(ek -e fi + E k ))(kl - (y^p - a(e k Q - e fj, - E k ) ((5 v e h ) j\^ pe 
h 1 i det Q h h 

e,h=— 1,1 " 

Y /^/«(e kp-efj,- E k ){kl - (uj 2 ) e h ) + a(e k - e \i - E k )r] e h jj^ e ^ e 

e,h=-l,l det Q h 

where: 

Ml = (E k + e /i) 2 + C e xr , (42) 
(wsOJ = (£ t + etf + Ja e , (43) 
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and: 



01 = Pft + aPift, (44) 

Xl = a~Pr h , (45) 

v\ = Pa + a 2 ^ e } (46) 

rf h = il-dr h e - (47) 

This propagator could be diagonalized but the diagonalization depends on 
k Q and and is quite complicated. It is not necessary to diagonalize the 
complete propagator for the purposes of this paper. 

Assuming that only the 75 condensate contributes, this corresponds to 
the relations: 

e h = xi = m = ehtx, ^ = o, v h = p& v e h = e h , m 

and it can be shown that: 

G + = V fco + e(E k - e /z) 

Y> fc + e(E fc - e /i) 

^^-(^-e^-ta) 2 7o ' 



which agrees with the quark propagator derived in [2^] using this ansatz. 



While this ansatz is a good first approximation, it will be shown that the gap 
equation will not close under this ansatz and the general quark propagator 



is the more complex form shown in fl4"I|). 



The analysis up until this point generalizes to perturbation theory. In 
what follows the analysis will be restricted to NJL models where things sim- 
plify considerably. With foresight define: 

n = fr-As, (50) 
a 

C h = h(eA 1 -aA 3 ), (51) 

where Ax is the usual CFL condensate and A 3 is a new condensate and the 
gap matrix is now given by: 

A = A l75 + A 3 7 57o . (52) 
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In this case: 



det n + = det = (jfeg - (£+) 2 ) (*; 2 - (e") 2 ) , (53) 
where the poles of the quasiparticle propagator are given by: 



(e ± ) 2 = E 2 +/i 2 + A 2 + A 2 ±V(2^/i + 2«A 1 A 3 )^ + 4/5 2 A 2 (E 2 + A 2 ) . (54) 

Note that e ± are distinct from the e ± defined in flSTp and are used exclusively 
in what follows. These dispersion relations can be shown to agree with those 
stated in [ETT . 



The poles of the propagator give the dispersion relations for quasipar- 
ticles, (-£~), quasi-antiparticles (s + ), quasiparticle holes (e~) and quasi- 
antiparticles holes (e + ). The physical manifestation of the condensates is a 
gap in the quasiparticle spectrum between the maximum of the quasiparticle 
branch and the minimum of the quasiparticle hole branch: 

V |mm) ( £ | max) 2£ \min- ('-*'-') 

This means that exciting a quasiparticle requires a minimum energy, (p. 
In the massless limit A3 = 0, a = 0, Ek — \k\ and therefore the reduce 



to e T defined in ||27|| . In this case the gap is simply (p = 2A% at k = //. 

In the massive case, A 3 alters the dispersion relations. The minimum of 
the quasiparticle hole spectrum occurs at: 



and leads to a gap: 



, 2 1 ,2 , A2 gm^Aj. + (A 2 - rrf) A 3 
^ +A 1 + A 3 A 2 + ^ 2 , (5b) 



^ = 2 /iAl - mA3 . (57) 



yU 2 + A 



2 



For A 3 = the minimum occurs at E k = fi and leads to a gap of <p = 2A X . 
When A 3 = is taken into account the minimum is shifted to a slightly 
lower value and the gap increases or decreases depending on the sign of A 3 . 
Further discussion of the dispersion relations and the gap will be given in 
subsequent sections in terms of specific solutions. 
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In the NJL model the full propagator is; 

G + (k) = (58) 
^ ae{k - fi + e E k ))(k% - MX) - e(k - /i - e E fc )/3 3 AiA 3 

ejklA (*8 ~ ( £+ ) 2 ) (*8 - (^) 2 ) 

P/a(e k -en- E k )(k 2 - {u 2 )%) + {k - fi - e £ fe )(Ai - eA 3 /a)A 3 ^ e 



(fco 2 - (^ + ) 2 ) {kl - (^) 2 ) 



where: 



Ml = (^ fc + e/i) 2 + (A 1 + eaA 3 )(A 1 + e A 3 /a), (59) 
(cu 2 )£ = (E fc + e/i) 2 + (A 1 + eaA 3 )A 1 . (60) 

The quasiparticle propagator can now be used in the mean field gap 



equation p7 1 



m = -ig 2 1 E r»D% B (k - q )G ( q )A( q )G+( q )r*, (6i) 

J A,B=1..N$-1 

where D^ B (k — q) is the gluon propagator, is the interaction vertex, 
and(see Appendix A): 

T = Cr T C _1 . (62) 

In the following sections results are presented for two different types of 
four fermion interactions. The next section presents results for the scalar 
NJL model. Section V presents results for an NJL model which has the color 
structure of single gluon exchange. These models give identical results if A 3 
is neglected but differ quite significantly if it is included. 

One of the main results of this section is the general quasiparticle propa- 
gator (|¥ID which could be used in a perturbation theory analysis of the gap 
equation. The second main result of this section is the quasiparticle prop- 
agator for an NJL model which is used in the rest of this paper. The 
poles of the quasiparticle propagator ([54]) are particularly important for this 
analysis. Finally the definition of the gap in terms of the condensates fl57|) is 
another result which will be important in the rest of this paper. 



13 



IV Scalar NJL Model 

The simplest possible four fermion interaction, motivated by the effective 
Lagrangian approach, is to take the interaction vertex to be = % and: 

9 2 DTb - j^jG 5^5 AB , (63) 
giving the matrix gap equation: 

A(k) = 8 i tt 2 Gj ^G (q)A( q )G + (q). (64) 

Acting on both sides of this equation with the operators 75, 7570 and 
tracing over the spinor indices one can obtain the gap equations for Ai and 
A3 respectively. The color-flavor structure of the gap matrix can almost be 
ignored because the propagators do not have any color flavor structure and 
therefore the same factor will occur on both sides of the equation. If one 
is working in the degenerate three flavor case there is a slight complication 
to even the zero mass gap equation because of the existence of the A 66 
condensate. This complication is ignored in this paper so this analysis is an 
approximation in the three flavor case. 

In the two flavor case there is only one color-flavor structure to be con- 
cerned about and this analysis is exact. This applies even to the case of the 
2SC color superconducting phase in the physical 3 flavor case as condensates 
involving the third color and flavor will simply decouple from the condensates 
of interest. 

The coupled gap equations are: 

J (2tt) 4 {ql - (e+)2) {ql - ( e -)2) 1 ' 



d A q 
(2?r 

2E q ^A 1 + A 3 (g 2 + A 2 - ^) + A 3 (/? 2 - a 2 )(£ 2 - A 2 ) 



A3 = -S^Gtt 2 /— ^ (66) 



where terms linear in q have been dropped since they will cancel out on 
integration over q from —00 to 00. The dependence of the condensate on 
momentum has been dropped since the right hand side of the gap equations 
are independent of k. 
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Acting on flB3| ) with the operators 70 and the identity and taking the 
trace leads to vanishing of the right hand side which is consistent with the 
assumption that A4 and Ag are zero. Acting with any of the other operators 
involving 7 ■ k will lead to a term involving q ■ k which will vanish by symmetry 
under the angular integration. 

Evaluation of these equations can be facilitated by the analytic contin- 
uation q — ► — zg 4 . The g 4 integration is then done by contour integration 
closing the contour in the upper half plane and picking up the poles at ie + 
and ie~ . The angular integrals can be done trivially giving: 



A 1= G 



dqq — + — 

V e + e~ 



4(m fi + AiA 3 )A 3 4(m fi + AiA 3 )A 2 



((.+) 2 - (e-r)s+ ((e+)> - (e-)2) e - 



(67) 



4 m /i Ai 




4m/i A] 



((.+) 2 - 



(68) 



If one assumes A 3 = the first equation reduces to the gap equation solved 



in [Q. The second equation does not vanish exactly under this assumption. 
If one takes m = and A 3 = the second equation is trivially satisfied and 
the first equation becomes the gap equation for massless quarks. 

The range of integration for q is not infinite since the NJL model is a 
four-fermion interaction model and must have an UV cutoff, A, which is left 
arbitrary except for the restriction that it be greater than fi. 

An approximate solution of the gap equation for m = is given by|J: 



Af } = Ai(m = 0) w 2 J A 2 -/i 2 exp 



A 2 - 
2/i 2 



exp 



2G/i 2 



(69) 



Expanding the integrands of (|67|) and (68) in m one can obtain approxi- 
mate solutions for the condensates: 



Ai(m)«Ai 



(o) 




C 1 - C(A 2 - fi 2 
[ 4G(i 4 
G(A 2 - fi 2 ) 
AG fi 4 



2 (A 2 - fi 2 
1 



+ A 



2(A 2 - fi 2 ) 



+ A 



1 - G(A 2 + fi 2 
2G/i 4 

l-G(A 2 + fi 2 ) ] 2 
2Gfi 4 J 
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3(G(A 2 - /i 2 ) - 1) 1 



16Gfi 6 4(A 2 - fi 2 ) 2 



+ 




A 3 (m) 




A — Ai(m) 



(70) 



(71) 



where: 



.4 



1 - 



3A 2 - /i 2 



(72) 



3A 2 -7/i 2 + 3/i 2 log (A^ 0) ) /4(A 2 -/i 2 ) 



In order to verify the accuracy of the approximate solutions the exact gap 
equations were also solved numerically using Mathematica © for comparison. 
These approximations are somewhat complicated due to their generality and 
range of applicability so some general comments are in order. 

The solution for A3 is linear in m//i for small values of m with slope A 
ranging from A — 1 for A^ = to A — for A^ = f a/A 2 — fj,' 2 . 

The solution for Ai(m) is an expansion in ^ to order 0(m 4 / /i 4 ). Factors 
of 1/Gfi 2 in some terms are actually enhancement factors as can be seen by 
taking typical values of G = 6 GeV _2 /7r 2 and [i = 500 MeV giving a value of 
1/Gn 2 w 6.5. Expanding only to 0(m 2 /fi 2 ): 



order of magnitude. This means that neglecting the effect of A3 by setting 
A = will have a non-trivial effect on the solution for Ai(m). The order 
(9(m 4 //i 4 ) terms are less instructive and are only included in order that the 
approximation is accurate for a wider range of m, //, A and G. 

The accuracy of this approximation for a range of four fermion couplings, 
G, is illustrated in Figure (jH). Examining the limiting case G — ► which 
corresponds to A^ — ■> 0, it can be seen that effect of the quark mass be- 
comes increasingly important. Conversely, as A^ increases the effect of the 
quark mass becomes less significant. Results obtained by varying A or fi 
have exactly the same behavior and also agree well with the exact numerical 
results. 
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Ai(m)/Ai(m=0) 




20 40 60 80 100 120 140 160 
m [MeV] 

Figure 1: Ai(m) for \l = 500 MeV and A = 2/i in the scalar NJL model for 
(from bottom to top) G7r 2 [GeV~ 2 ] = 3.5,4.875,6. Analytic approximations 
and exact numerical results are shown for comparison. 

The solutions for G = 3.5GeV~ 2 /7r 2 and G = 6GeV~ 2 /n 2 which corre- 
spond to A5 0) = 10.14 MeV and AS 0) = 103.27 MeV are shown in Figures 
d) and (|) with // = 500 MeV. Also shown in Fi gure (H) are the solutions 
where A3 is neglected for comparison. The effect on Ai(ra) of including A3 
is an increase in the mass dependence compared to the solution neglecting 
A3. The effect is more pronounced for smaller values of Aj°\ but is still sig- 
nificant at A^ « 100 MeV. Notice that the relative size of A3 with respect 
to Ai is very similar in both cases. 

As discussed in the previous section, inclusion of A 3 alters the dispersion 
relations and the gap, <p. The gap between the quasiparticle hole branch and 
quasiparticle branch in this case is not simply <p = 2Aj(m) but is given by 
(|5"7|). The value of the gap is plotted for the same choices of G presented 
above in Figure with and without A3. The effect on the gap is larger 
than the effect on Ai alone because A 3 affects (p directly in (|57|) as well as 
indirectly through Ai. The dispersion relation for the quasiparticle holes is 
not appreciably altered. 
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Ai (m) 
Ai (m=Q 




25 50 75 100 125 150 25 50 75 100 125 150 

mlMeV] m^leV] 



Figure 2: Ai(m) for [i = 500 MeV and A = 2fi in the scalar NJL model. 
The graph on the left is for G = 3.5 GeV~ 2 /7r 2 and A^ 0) = 10.14 MeV. The 
graph on the right is for G = 6 GeV~ 2 /7r 2 and aS 0) = 103.27 MeV. Solutions 
are shown neglecting and including A 3 . 
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0.2 
0.175 
0.15 
0.125 

A 3 (m)/Ai(m=0) 0.1 
0.075 
0.05 
0.025 

20 40 60 80 100 120 140 160 
m [MeV] 

Figure 3: A 3 (m) for /i = 500 MeV and A = 2/x in the scalar NJL model. 
Four fermion coupling constants G = 3.5 GeV~ 2 /7r 2 and G = 6 GeV~ 2 /7r 2 
were used and give solutions with AS 0) = 10.14 MeV and = 103.27 
MeV. Analytic approximations and exact numerical results are shown for 
both cases. 



A 1 (m=0)=10.14MeV 
A!(m=0)=103.27 MeV 

0.69 m/|iAi(m) 
m/2|i Ai(m) 
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25 50 75 100 125 150 '0 25 50 75 100 125 150 



m [MeV] m [MeV] 

Figure 4: Plot of tp/tfo for the scalar NJL model including and neglecting A 3 
for different values of Aj°\ The graph on the left is for G = 3.5 GeV~ 2 /7r 2 
and aS 0) = 10.14 MeV. The graph on the right is for G = 6 GeV^/n 2 and 
AS 0) = 103.27 MeV. 
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V NJL Model motivated by Single Gluon Ex- 
change 

The NJL interaction which has the structure of single gluon exchange has 
the interaction vertices: 

ri = 7,f, r£ = -7m(¥) t , (74) 

and: 

g 2 D^ B -> 3tt 2 G^^ ab , (75) 

giving the gap equation: 

d A q 



= -2 i G?/ ^ 7 , G (g)A(g)G+(g) 7", (76) 
using the equations: 

g 

(A )^A^Ajy = — assuming \ B is antisymmetric. (77) 



Multiplying each side of equation ( [75p by 75 and 7570, tracing over the Dirac 
indices and using the cyclicity of the trace and the relations: 

7*tV = -47 5 , (78) 

7m7 5 7V = -7 5 7m7V = 2 7 5 7°, (79) 

and then applying all the other machinery of the last section, the coupled 
gap equations are obtained: 



A 1 = G j ' dqq 2 ^ 



Ai , A x 4(to h + AiA 3 ) A 3 4(m y, + A X A 3 ) A 3 



A, 



e+ e- ((e+) 2 - (.-) 2 ).+ ((.+) 2 - (.-) 2 ).- J ' 

(80) 

G r 2 ( A 3 A 3 4 m /j, Ax 4m/iAi 



2i ^ \ e+ e- ((e+) 2 - (e-) 2 )e+ ((e+) 2 - (e-) 2 y 



4(g 2 + A 2 )A 3 4 (g 2 + A 2 ) A 



3 



((£+)' - (.-) 2 ) £ + ((.+) 2 - (.-) 2 ).- 



(81) 
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Notice that they are almost the same as the equations in the last section 
except for an overall factor of — ~ in the second equation. Therefore the only 
change to the approximate solutions of the last section is: 



1 3A 2 - V 
A = -- + 



z 3A 2 - 10/i 2 + 3/i 2 log (aS 0) ) /4(A 2 -/i 2 ) 
For clarity the solution to order 0(m 2 / /j?) is repeated here: 
A 1 (m) = Af ) 



(82) 



^ 2 f l-g(A 2 -^ 2 ) _ 1 + A \ - G(A 2 + /i 2 ) 



4G> 4 2(A 2 -/x 2 ) 2GV 

(83) 

Tfi 

A 3 (m) = A — Ai(m). (84) 
A 4 

The accuracy of this approximation for a range of these parameters is illus- 
trated in Figures ©-© by comparison with exact numerical solutions. 

The main difference from the previous section is in the sign of A. A 3 (m) 
is again linear in m//j for small m but with a slope of A « — | for small 

values of A^ and approaching as A^ increases. 

The first and third terms in the 0(m 2 /n 2 ) expansion of Ai(m) are of the 
same order of magnitude but in this case their signs differ. The resulting effect 
of the inclusion of A 3 is to reduce the relative effect of the mass dependence. 

The comments on the mass dependence as a function of the parameters 
jj, and G are unchanged from the previous section. For small values of A^ 
the effect is important and is less significant for increasing Af } . This effect 
can be seen in Figure ([5]). 

In the weak coupling limit, G — > 0, equations (|83| ) and (|34D are very 
similar to equations (37) and (34) of 



The mass dependence of Ai as a function of A is different in this model 
than in the last section. In this case the last two terms of equation ( j83|) 
decrease as A approaches /z and although the first term still dominates, the 
relative mass dependence decreases with decreasing A. The effect is that the 
mass dependence is less significant for smaller A^° as can be seen in Figure 

(§)■ 

The solutions for G = 3.5 GeV 2 /7r 2 and G = 6 GeV 2 /7r 2 which cor- 
respond to Ai(m = 0) = 10.14 MeV and A x (m = 0) = 103.27 MeV are 
shown in Figures (0) and @ with fi = 500 MeV. The solutions where A 3 is 
neglected are also shown in Figure (^) for comparison. The effect on Ax(m) 
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of including A 3 is a reduction in the relative effect compared to the solution 
neglecting A 3 . The effect is more pronounced for smaller values of Aj°\ but 
is still nonneglegible at A^ pa 100 MeV. Notice that the relative size of A 3 
with respect to Ai is very similar in both cases. 

Inclusion of A3 alters the dispersion relations and the gap. The gap 
between the quasiparticle hole branch and quasiparticle branch is not simply 
<p = 2Ai, but is again given by fl5"T|). The gap depends on A 3 both explicitly 
and implicitly through its effect on A^ The value of the gap is plotted for 
the same choices of G presented above in Figure (|S]) with and without A 3 . 
The effect on the gap is larger than the effect on Ai alone. The zeroth order 
approximation, Ai(m = 0), is a better approximation to the gap than the 
solution where A 3 is neglected. The dispersion relation for the quasiparticle 
holes is not appreciably altered. 

Extending this analysis out to quark mass of m — 150 MeV is partly for 
illustrative purposes but is also valuable for two other reasons. 

First the analysis of |22[ shows that in a coupled analysis of the super- 
conducting (diquark) condensate and the axial condensate, constituent quark 
masses of the order of 100 MeV for the light quarks are possible at fi pa 400 
MeV. For constituent quark masses of this order the results of our analysis 
show that the presence of a new condensate, A 3 which they neglected in their 
ansatz will be relevant. Their ansatz is a reasonable first approximation as 
the new condensate and its effects on Ai are not large. However, the ef- 
fect on the gap ip itself is significant even for gaps of the order of 100 MeV. 
The full analysis requires the more general ansatz and the use of the general 
quasiparticle propagator. The effects of chiral symmetry breaking have not 



been taken into account in this paper as in E4 21] but such an extension 



would not be difficult. The authors of |^TJ also found significant values for 
the constituent quark masses for fi pa 400 MeV. 

Second the methods used in this analysis can be extended to the physical 
case where the strange quark mass is of the order of 150 MeV. It is instructive 
therefore to carry out this analysis as a precursor to the physical case. 

The two NJL models analyzed in this paper lead to quite different results 
when the full calculation is performed. In Figure (pi]) it is shown that the 
full calculation depends significantly on the type of NJL model used. For 
lower values of A^ the effect on the gap is almost a factor of 2 and is still 
significant even for A^ of the order of 100 MeV. The second NJL model 
analyzed is perhaps more physically motivated than the second but there is 
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20 40 60 80 100 120 140 160 
m [MeV] 

Figure 5: Ai(m) for (j, = 500 MeV and A = 2/i in the NJL model 
with the color structure of single gluon exchange for (from bottom to top) 
G7r 2 [GcV~ 2 ] = 3.5,5.25,6.5. Analytic approximations and exact numerical 
results are shown for comparison. 
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Ai(m)/Ai(m=0) 
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■ Ai(m=0)=103.27MeV 




★ Ai(m=0)=48.20MeV 




♦ Ai(m=0)=28.60MeV 





20 40 



60 80 100 
m [MeV] 



120 140 160 



Figure 6: A^m) for /i = 500 MeV and G = Q GeV^ 2 / 7 ^ in the NJL model 
with the color structure of single gluon exchange for (from bottom to top) 
A = 2yU, 1.7/1, 1.5/1. Analytic approximations and exact numerical results are 
shown for comparison. 
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Figure 7: Ai(m) for fi = 500 MeV and A = 2/x in the NJL model with 
the color structure of single gluon exchange. The graph on the left is for 
G = 3.5 GeV- 2 /ir 2 and a! 0) = 10.14 MeV. The graph on the right is for 
G = 6 GeV- 2 /n 2 and AS 0) = 10.14 MeV. Solutions are shown neglecting and 
including A 3 . 
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Figure 8: A 3 (m) for fi = 500 MeV and A = 2/x in the NJL model with the 
color structure of single gluon exchange. Four fermion coupling constants 
G = 3.5 GeV~ 2 /7r 2 and G = 6 GeV~ 2 /7r 2 were used and give solutions with 
A^ 0) = 10.14 MeV and A^ 0) = 103.27 MeV. Analytic approximations and 
exact numerical results are shown for both cases. 




-0.27 m/\L Ai (m) 



27 



(p/(po 




25 50 75 100 125 150 
m [MeV] 



cpApo 



1 

0.975 
0.95 

0.925 
0.9 

0.875 
0.85 

0.825 



neglecting A 3 
including A3 



25 50 75 100 125 150 
m [MeV] 



Figure 9: Plot of <p/(po for the NJL model with color structure of single gluon 
exchange, including and neglecting A3 for different values of Af\ The graph 
on the left is for G = 3.5 GeV" 2 /^ 2 and A i 0) = 10 - 14 MeV - The g ra P h 011 
the right is for G = 6 GeV- 2 /?r 2 and A^ 0) = 103.27 MeV. 
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25 50 75 100 125 150 '0 25 50 75 100 125 150 



m [MeV] m [MeV] 

Figure 10: Comparison of the full gap calculations in the scalar NJL model 
and the single gluon exchange motivated NJL model. The graph on the left 
is for G = 3.5 GeV" 2 /^ 2 and A i° } = 10.14 MeV. The graph on the right is 
for G = 6 GeV- 2 /7i 2 and A^ 0) = 103.27 MeV. 
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is no reason to restrict the analysis to interactions of this form. The most 
general NJL model would involve both types of interactions. This analysis 
shows that there could be a significant dependence on the specific NJL model 
used. 



VI Conclusion 

In this paper the general quasiparticle propagator for the case of equally 
massive quarks is derived. The quasiparticle propagator is then specialized 
to the propagator in an NJL model. It is shown that there is exactly one 
new condensate, A 3 , in this model. 

The quasiparticle propagator was used to solve the gap equations for the 
CFL condensate, A 1; and the new condensate, A 3 in two NJL models as a 
function of the quark mass. Approximate analytic solutions were obtained 
in both cases valid over a range of values of the parameters m, A, /i and 
G in both models. The accuracy of these approximations was tested by 
numerically solving the exact gap equations. This represents a significant 
advance as the analytic approximations render numerical solution of the gap 
equations unnecessary for a significant range of parameter space. 

Results for the condensates, Ai and A3, and the gap, (p, in the scalar NJL 
model were presented as a function of the quark mass. The full solution for 
Ai and (p were compared to solutions obtained neglecting A 3 . This was done 
for different values of the four fermion coupling constant. The results show 
that the complete solution differs from the approximate solution especially 
at large quark mass and smaller Ai. 

Results for the condensates, Ai and A 3 , and the gap, <p, in the NJL model 
with the structure of single gluon exchange were presented as a function of 
the quark mass. The full solution for Ai and (p were compared to solutions 
obtained neglecting A3. This was again done for different values of the 
four fermion coupling constant. The results show that the complete solution 
differs from the approximate solution especially at large quark mass and 
smaller A^ 



The solutions obtained in this paper are relevant to the work of [^] where 
A 3 was neglected as a first approximation. The approach of this paper could 
be used to generalize their analysis to the complete case. 



The authors of pTJ showed that A 3 could be quite large in an NJL model 



with a Lorentz non-invariant interaction and have a significant effect on the 
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gap. In the analysis presented here it is shown that even for small values of 
A3, neglecting this condensate can have a non-trivial effect on the gap. 

Comparison of the results for the gap in the two different NJL models 
analyzed in this paper show that the gap is strongly dependent on the type 
of NJL model used. 

Solving the color superconducting gap equations for the case of equal mass 
quarks is also valuable as a precursor to analyzing the physical case where 
the up and down quarks are essentially massless and the strange quark is 
massive. The methods used in this research combined with the methods of 



T5f can be combined to analyze the physical case which is the ultimate goal 



of this line of research. 
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A Nambu-Gorkov Formalism 

In the standard Nambu-Gorkov formalism one deals with the eight compo- 
nent spinors: 

* = ( ^ ) , *=Nc), (Al) 

where C is the charge conjugation operator. The action for interacting 
fermions in the mean field approximation can be written concisely as: 

= - J $ 5-1 ^ (A2) 

where: 

s ' 1 - ( |g aT [<£]""-. ) • <*>) 

Go are the free propagators of particles and charge-conjugate particles, A + 
is a matrix of diquark condensates and A - = 7oA + 7q. 
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The full Nambu-Gorkov propagator is: 

S ={-G Q A+G+ G- J' < A4 > 

where: 

G* = {[(ft]- 1 - XrcfiA*}- 1 . (A5) 
are the full quasiparticle and charge-conjugate quasiparticle propagators. 



The gap equation is a consistency equation p?fl: 
A + (k) = ~ig 2 j 7^4 E r£D% B (k-q)G^ q )A + ( q )G+( q )T*. (A6) 



A,B=1..N?-1 



which must be satisfied in order that the effects of interaction at the mean 
field level are consistently incorporated in the action ( |A2|) in terms of the 



diquark condensate matrix^, A + . D^ B {k — q) is the gluon propagator, is 
the interaction vertex, and: 



T = CT T C-\ (A7) 



B Dirac Basis 



The Bailin and Love motivated decomposition given in 27 



is: 



A = A175 + A 2 7 ■ &7075 + A37075 +A 4 + A57 ■ lcy + A 6 7 ■ k + A77 ■ k^ 5 + A 8 7 , 

(Bl) 

where: 

A175 + A 2 7 ■ kjolb, (B2) 

represents condensation of fermions with the same chirality in the even parity 
channel, 

A 4 + A57 ■ ho, (B3) 

represents condensation of fermions with the same chirality in the odd parity 
channel. 

A 3 7o75 + A 7 7 • Ary 5 , (B4) 

4 A + is also called the gap matrix since it elements determine the gaps in the quasipar- 
ticle spectrum. 



32 



represents condensation of fermions with opposite chirality in the even parity 
channel, 

A 6 7 • k + A 87o , (B5) 

represents condensation of fermions with opposite chirality in the odd parity 
channel. 
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